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Abstract 

We report the results of a precise measurement of the D s meson lifetime based 
on 1662 ± 56 fully reconstructed D s — ► (fin decays, from the charm hadroproduction 
experiment E791 at Fermilab. Using an unbinned maximum likelihood fit, we mea- 
sure the D s lifetime to be 0.518 ± 0.014 ± 0.007 ps. The ratio of the measured D s 
lifetime to the world average D° lifetime [1] is 1.25 ± 0.04. This result differs from 
unity by six standard deviations, indicating significantly different lifetimes for the 
D s and the D°. 

PACS: 13.25.Ft, 14.40.Lb 



Precise measurements of the lifetimes of the weakly decaying charm mesons are 
useful for understanding the contributions of various weak decay mechanisms. 
Despite the fact that they are all tied to the charm quark decay, the decays 
of the ground-state charm mesons can have different contributions from the 
four first-order processes (two spectator, W- annihilation, and W-exchange), 
and the lifetimes are, in fact, quite different [1]: 

t(D + ) : t(D°) : t(D s ) = 2.5 : 1 : 1.1 



The least well measured among these lifetimes is that of the D s , which is known 
only to about the 4% level [2]. It is clear that there is a very large difference 
between the lifetime of the D + and that of the D° or D s . However, the current 
measurements are only suggestive of a difference between the lifetimes of the 
D° and the D s : 

T |^| = 1.13 ± 0.04 ( 3<r difference from unity ). 
r(.D u ) 



A more precise measurement of the lifetime of the D s could more clearly 
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establish whether the lifetimes of the D° and D s are indeed different. 

The large difference between the D + and D° lifetimes might be explained by 
a large contribution from W-exchange in D° decays and/or large destructive 
interference between D + decay amplitudes. A difference between the D° and 
D s lifetimes could be due to phase space differences, color factors, W-exchange, 
W-annihilation, decay via r lepton modes, or relativistic effects [3,4]. The 
sensitivity of the D°-D s lifetime difference to W-annihilation is noted in the 
paper by Bigi et al. [4]. However, spectator decays of the D s have one more 
strange quark than spectator decays of the D°; hence, the Cabibbo-favored 
hadronic final states have an extra kaon and correspondingly smaller phase 
space. Table 1 compares the phase space available for 2, 3, 4, and 5-body 
decays of the two mesons. The table shows that phase space differences alone 
could account for the currently observed lifetime differences. 

In this letter, we report results of a high-statistics measurement of the D s 
lifetime based on data from the charm hadropro duct ion experiment E791 at 
Fermilab. (Unless otherwise specified, charge conjugate states are implicitly 
included in this letter.) This measurement establishes more certainly a dif- 
ference between the lifetimes of the D° and the D s . We use only the decay 
mode D s — > (fin since the (fi mass constraint gives a D s sample with a very good 
signal-to-noise ratio. The reconstruction efficiency as a function of lifetime was 
determined from data for the Cabibbo-favored (CF) decay D + — > K~n + n + 
and from the ratio of the efficiency for D s — > (fin to that for D + — > K~n + n + 
obtained from Monte Carlo (MC) simulations. With this technique, most un- 
certainties that arise from MC simulation should cancel. 

E791 [5] is a high statistics charm experiment which acquired data at Fermilab 
during the 1991-1992 fixed-target run. The experiment combined a fast data 
acquisition system with an open trigger. Over 2 x 10 10 events were collected 
with the Tagged Photon Spectrometer using a 500 GeV 7r~ beam. There were 
five target foils with 15 mm center-to-center separations: one 0.5 mm thick 
platinum foil followed by four 1.6 mm thick diamond foils. The spectrometer 
included 23 planes of silicon microstrip detectors (6 upstream and 17 down- 
stream of the target), 2 dipole magnets, 10 planes of proportional wire cham- 
bers (8 upstream and 2 downstream of the target), 35 drift chamber planes, 2 
multi-cell threshold Cerenkov counters that provided ir/K separation in the 
6-60 GeV/c momentum range [6], electromagnetic and hadronic calorimeters, 
and a muon detector. 

All criteria used to select candidate D s — > (fin decays, with the exception of 
Cerenkov identification requirements and the (fi mass cut, were chosen to max- 
imize S/y/S + B, where S is the number of D + — * K~n + n + signal events in 
data, scaled to the level expected for D s — > (fin, and B is the number of back- 
ground events from appropriate sideband regions of the (fin mass distribution. 
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To reconstruct the D s — > (fin and D + — > i^~7r + 7r + candidates, three-prong 
decay vertices with charge of ±1 were selected. All decay tracks were required 
to travel through at least one magnet and be of good quality. Decay vertices 
were required to be located outside the target foils, and the significance of 
spatial separation from the primary vertex in the beam direction, Az/<7az, 
where <ta z is the error on the separation, was required to be at least 11. The 
component of the D momentum perpendicular to the line joining the primary 
and secondary vertices was required to be less than 0.35 GeV/ c. The transverse 
impact parameter of the D momentum with respect to the primary vertex was 
required to be less than 55 /im. Decay tracks were required to pass closer to 
the secondary vertex than to the primary vertex, and the sum of the squares 
of their momenta perpendicular to the D direction was required to be larger 
than 0.3 (GeV/c) 2 . The candidates in the decay mode (fin with (fi — > K + K~ 
were also required to have M(K + K~) within ±10 MeV/c 2 of the (fi mass. 

For the decay D + — > K~n + n + , the kaon was identified on the basis of charge 
alone and no Cerenkov identification cuts were needed. For the two same-sign 
tracks in the decay D s — > (fin, the track with the highest Cerenkov probability 
was assumed to be the kaon and the other track the pion. No further Cerenkov 
identification cuts were applied, given the (fi mass selection criteria. 

The longitudinal and transverse position resolutions for the primary vertex 
were 350 /xm and 6 //m, respectively. The mean momentum of selected D + 
and D s mesons was 70 GeV/c. For the decay vertices of these mesons, the 
transverse resolution was about 9 /iin, nearly independent of the momentum, 
and the longitudinal resolution was about 360 /im at 70 GeV/c and worsened 
by 30 /im for every 10 GeV/c increase in momentum. 

To eliminate any possible background due to reflections from D + — > K~n + n + 
with a pion misidentified as a kaon, which appear under and near the D s signal 
in the (fin mass plot, we excluded all (fin candidate events with M(K~n + n + ) 
within ±30 MeV/c 2 of the D + mass. The mass distributions for all (fin can- 
didates that survived our selection criteria are shown in Fig. 1. The higher 
mass peak corresponds to D s — > (fin while the lower peak corresponds to the 
singly Cabibbo-suppressed (SCS) decay D + — > (fin + . The (fin mass distribu- 
tion before D + — > K~n + n + background subtraction is shown in Fig. 1(a), 
with candidates consistent with D + — > K~n + n + background shown as the 
hatched region. The (fin mass distribution after the reflection subtraction is 
shown in Fig. 1(b). The hatched region in Fig. 1(a) appears to be due to a 
small D s signal and a linear background that turns on around 1.95 GeV/c 2 . 
The background in Fig. 1(b) is therefore modeled as a piecewise linear func- 
tion with a discontinuity fixed at 1.95 GeV/c 2 . We have studied the sensitivity 
of our measured lifetime to the position of this discontinuity and have found 
very little effect. We have included a small systematic error in our total error 
to account for the change in lifetime when the location of this discontinuity is 
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varied from 1.95 to 1.99 GeV/c 2 . 



To extract the D s lifetime we performed a simultaneous unbinned maximum- 
likelihood fit to the distribution of mass and reduced proper decay time (t R ) 
of all (fin candidates that pass the selection criteria. The t R is defined as 

t R = (L — L roin )/c/3 7 , 



where L is the distance between the secondary and primary vertices, and 
L min is the minimum L value allowed by the cut on the separation between 
the secondary and primary vertices for this event. We fit events in the ranges 
1.75 < M(0tt) < 2.05 GeV/c 2 and 0.0 < t R < 4.0 ps. From the fit we also 
measure the lifetime of the SCS decay D + — > <pn + . 

The overall likelihood function is 

p (,-^pred ^obs ) l^-^pred 

c=- — 7— = — riA, 

v/27T Npred j 



with N pred = N D + + Nd s + Nbg> where N D +, Nd s , and Nbg are the fitted 
number of D + — > (fin + , D s — > (fin, and (fin background events, respectively. 

The likelihood for each candidate is 

d = SD+irrn, tf ) + S Ds (mi, t R ) + B(m t , t R ) 



where 



S D+ {m h tf) = F D+ x -=L e ~{™ D+ ~^) 2 /^ D+ '' 

xR D+ -f D+ {t R )-e' l ?^\ 



S Ds ( mi ,t R ) = F Ds x 



e ~{m Ds -mi) /2a Ds 

2na Ds 

* Rd s ■ f Ds (t?) ■ e- t?/TD % 



B(m t ,t R ) = F BG x R BG ■ ( e -^ /TB ^ + Ce'^'^) 

A 1 + Siirii ( for rrii < m ) 

< 



A 2 + S 2 rrii ( for m, > m ) 
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The functions Sb+(ith, tf), Sd s (ith, tf), and B(rrii, tf) are the likelihood func- 
tions for D + — > 07r + , D s 07r, and <f)n background, respectively. The coeffi- 
cients Fe>+, Fd 3 , and Fbg are the ratios of the D + , D s , and the background, 
to the total number of events: F D + = N D +/N pred , F Da = N Ds /N pred , and 
Fbg = NBG/N pre d- rrii and tf are the mass and the reduced proper decay 
time of the (pir candidates. &£>+ and gd s are the D + and D s mass resolutions. 
td+ and Td s are the D + and D s lifetimes. The background in the mass dis- 
tribution is described as two linear functions with slopes Si and S2 which are 
determined from the fit, and m is the position of the discontinuity in the 
background mass shape, which is fixed. The terms A\ and A 2 are obtained 
from the normalization and the size of the discontinuity, which is determined 
from the fit. The coefficients Rp+, Rd s , and Rbg are the normalizations of 
the fit functions over t R in the limits of the fit. 

The functions f D + (tf ) and f Ds (tf) are the reconstruction efficiency functions 
of t R , for D + — > 07r + and D s — > 07r. These efficiency functions are determined 
from the CF D + — > K~7t + tt + decay in data, and from the ratios of MC 
efficiencies for the SCS D + — > 07T + decay and the D s — > 07r decay, separately, 
to that for D + — > K~-n +r n + decay, as follows: 



f D+ (tf) = e D D ^ K ^(tf)* 
fD a (tf) = ^t K ^t R )x 



MC ( f R\ 

Tmc 77r~\i 

MC UR\ 

e D s ^<t>TT\ l i ) 

f MC I fR \ ■ 



We determine the efficiency for D + — > K^-k + -k + from data by comparing the 
background subtracted t R distribution for D + — > K~h + 'k + to an exponential 
function with the world average D + lifetime [1], r(D + ) = 1.057 ± 0.015 ps. 
The uncertainty on the D + lifetime contributes to our systematic error. 

We prefer to determine the functions fo(tf), fo+(tf) and fD a (tf), from data 
and the ratios of efficiencies from MC to reduce our dependency on MC. Most 
uncertainties that arise from MC simulation should cancel in the ratios and, 
as a systematic check, we also measure the D s lifetime using the functions 
fo+(tf) and fo s (tf) derived only from MC. We use the same cuts for the 
07r and K~7r + n + data samples in order to minimize any bias in the efficiency 
ratios due to the selection criteria. 

The background in the t R distribution is parameterized as a sum of two ex- 
ponential functions with lifetimes of tbg\ and t B g 2 > and a relative coefficient 
C. The parameters were determined from the t R distribution for the D s side- 
bands. As a systematic check they were also allowed to float in the unbinned 
maximum-likelihood fit, with no significant effect on the measured D s lifetime; 
a small systematic error is included in our total error to account for this effect. 
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The projections of the data and the best fit function from the unbinned 
maximum-likelihood fit on the M(07r) and t R distributions for all events that 
pass the selection criteria are shown in Figs. 2 and 3(a), respectively. For 
the subset of events within ±25 MeV/c 2 of the D s mass, the projections of 
these fit results on the t R distribution, scaled for the new mass range, are 
shown in Fig. 3(b). This illustrates that our fit describes the data well for 
events in a narrow window around the D s signal where the background is 
small. The fit yields 1662 ± 56 D s — > 4>n signal events with a D s lifetime 
T (D S ) = 0.518 ± 0.014 ps, and 997 ± 39 D + — > 0tt+ signal events with a D+ 
lifetime r(D + ) = 1.065 ±0.048 ps. The D + lifetime is in good agreement with 
current D + lifetime measurements, indicating that the MC models the relative 
efficiency of 07T to K~7t + tt + well. The projections from the fit on the t R distri- 
bution for D + — > 07r + and D s — > 07r are not pure exponentials because of the 
nonuniform detector efficiencies, which were parameterized by the functions 
fo+(t R ) and fD s {tf ) as described above. 

Since ratios of MC efficiencies are used in the fit, most systematic errors cancel 
in the final D s lifetime measurement. However, some uncertainties remain. The 
sources of these are summarized in Table 2, and quantitative estimates of the 
errors are listed. In general, the estimates come from re-fitting after varying 
data selection criteria, parameters, and mass or t R regions used. The total 
systematic error comes from adding the contributions in quadrature. 

Other checks were performed such as fitting to only a portion of the t R range, 
which would check our efficiency functions across the t R range considered, and 
measuring the lifetime of the and D~ separately. Variations in results were 
found to be small compared to the statistical error. Uncertainties due to target 
absorption and scattering of charm secondaries are expected to be small due 
to the thinness of the target foils and the requirement that decay vertices be 
located outside the target foils. In addition, since we use the D + — > K~7t + tt + 
data to correct for efficiency, any such effects should be accounted for. 

Finally, as a test of our fitting procedure, we fitted our K~it + it + data to 
measure the D + lifetime using an unbinned maximum-likelihood fit similar to 
that used to fit the 07r data. The efficiency function for the D + — > K~7t + tt + 
was determined only from MC. The fit yielded 62651 ± 337 D + signal events, 
shown in Ref. [7], with a measured D + lifetime that agrees with the world 
average value [1] to better than the 1.5% error on the world average value. 

In conclusion, we have made a new precise measurement of the D s lifetime us- 
ing 1662 ±56 fully reconstructed D s — > (pre decays. The D s lifetime is measured 
using an unbinned maximum-likelihood fit to be 0.518 ±0.014 ±0.007 ps. This 
value is about two standard deviations higher than the present world average 
of 0.467 ±0.017 ps [1]. Using our result and the world average D° lifetime [1], 
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we find the ratio of our D s lifetime to the D° lifetime to be 

t(D s . 



t(L>°) 



1.25 ± 0.04 (6<T difference from unity 



showing significantly different lifetimes for the D s and D°. This result may 
be used to constrain the contributions of various decay mechanisms to charm 
decay and further refine our quantitative understanding of the hierarchy of 
charm particle lifetimes. 
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Fig. 1. The <pir mass distributions: (a) before excluding candidates with 
) within ±30 MeV/c 2 of the D + mass, which are shown as the hatched 
histogram; (b) after excluding candidates with M(K~ir + 7r + ) within ±30 MeV/c 2 
of the D + mass, as described in the text. 
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Fig. 2. The projection from the unbinned maximum-likelihood fit on the M(c/)tt) 
distribution, with the discontinuity in the linear background fixed at 1.95 GeV/c 2 . 
The two peaks are the SCS D + signal with a fitted yield of 997 ± 39 events, and 
the D s signal with a fitted yield of 1662 ± 56 events, as described in the text. 
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Fig. 3. The projection of the data and the best fit function from the unbinned 
maximum-likelihood fit, on the reduced proper decay time distribution, (a) All 
events that pass the selection criteria (histogram). Also shown are the total fit 
function (solid curve) , and the contribution to the fit from D s — > (J)tt (dashed curve) , 
D + — > 4>ir + (dashed-dotted curve), and <frir background (dotted curve), (b) Events 
within ±25 MeV/c 2 of the D s mass (histogram). Also shown are the total fit function 
(solid curve), and the contribution to the fit from D s — > (f>ir (dashed curve) and from 
(pTr background (dotted curve), as described in the text. 
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Table 1 

Summary of the phase space (PS) factors, in arbitrary units, and the phase space 
ratios, for 2, 3, 4, and 5-body decay modes of the D° and D s . The phase space 
factors are calculated for the final states shown, with no intermediate resonances. 



D° 




D s 






Decay Mode 


PS 


Decay Mode 


PS 


PS(D°)/PS(D S ) 


K-1T+ 


1.45 


K+K° 


1.36 


1.07 


K-7T+7T 


1.70 


K+K--K+ 


1.32 


1.29 




0.537 


K°K-TT + TT + 


0.247 


2.17 


-w + + - - 


0.0548 




0.0127 


4.31 



Table 2 

Summary of the systematic uncertainties assigned to the D s lifetime, as described 
in the text. 



Source of Uncertainty Systematic Error 



D + lifetime uncertainty [1] 


0.004 


ps 


Size of D + — ► K~tt + tt + reflection 


0.003 


ps 


Data selection criteria 


0.003 


ps 


Likelihood efficiency functions, fo{tf) 


0.003 


ps 


t R parameterization of background 


0.002 


ps 


Position of step in background function 


0.002 


ps 


Monte Carlo production model 


0.002 


ps 


Total systematic error 


0.007 


ps 
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